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Abstract 



The Dirac equation is exactly solved for a pseudoscalar linear plus Coulomb- 
like potential in a two-dimensional world. This sort of potential gives rise to 
an effective quadratic plus inversely quadratic potential in a Sturm-Liouville 
problem, regardless the sign of the parameter of the linear potential, in sharp 
contrast with the Schrodinger case. The generalized Dirac oscillator already 
analyzed in a previous work is obtained as a particular case. 



The Schrodinger equation with a quadratic plus inversely quadratic po- 
tential, known as singular oscillator, is an exactly solvable problem pP-jH] 
which works for constructing solvable models of N interacting bodies [7j- 
[Hj as well as a basis for perturbative expansions and variational analyses 
for spiked harmonic oscillators 0-^1] • Generalizations for finite-difference 
relativistic quantum mechanics ^H] as well as for time-dependent param- 
eters in the nonrelativistic version have also been considered |ini-CZ|- In 
the present paper we approach the time-independent Dirac equation with a 
conserving-parity pseudoscalar potential given by a linear plus a Coulomb- 
like potential. This sort of potential gives rise to an effective quadratic plus 
inversely quadratic potential in a Sturm-Liouville problem, regardless the 
sign of the parameter of the linear potential. Therefore, this sort of problem 
is exactly solvable in the Dirac equation. The generalized Dirac oscillator 
[THj is obtained as a particular case. In addition to its importance as a new 
solution of the Dirac equation this problem might be relevant to studies of 
confinement of neutral fermions by a linear plus inversely linear electric field. 

The two-dimensional Dirac equation can be obtained from the four-dimen- 
sional one with the mixture of spherically symmetric scalar, vector and 
anomalous magnetic-like (tensor) interactions. If we limit the fermion to 
move in the x-direction {py = Pz = 0) the four- dimensional Dirac equation 
decomposes into two equivalent two-dimensional equations with 2-component 
spinors and 2x2 matrices Then, there results that the scalar and vector 
interactions preserve their Lorentz structures whereas the anomalous mag- 
netic interaction turns out to be a pseudoscalar interaction. Furthermore, 
in the 1+1 world there is no angular momentum so that the spin is ab- 
sent. Therefore, the 1+1 dimensional Dirac equation allow us to explore 
the physical consequences of the negative-energy states in a mathematically 
simpler and more physically transparent way. The confinement of fermions 
by a pure conserving-parity pseudoscalar double-step potential PD] and their 
scattering by a pure nonconserving-parity pseudoscalar step potential 
have already been analyzed in the literature providing the opportunity to 
find some quite interesting results. Indeed, the two-dimensional version of 
the anomalous magnetic-like interaction linear in the radial coordinate, chris- 
tened by Moshinsky and Szczepaniak j^] as Dirac oscillator, has also received 
attention. Nogami and Toyama Toyama et al. |21] and Toyama and 
Nogami 25^ studied the behaviour of wave packets under the influence of 
such a conserving-parity potential, Szmytkowski and Gruchowski [2^1 proved 
the completeness of the eigenfunctions and Pacheco et al. j27j studied some 
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thermodynamics properties. More recently de Castro [TH] analyzed the pos- 
sibility of existence of confinement of fermions under the influence of pseu- 
doscalar power-law potentials, including the case of potentials unbounded 
from below, and discussed in some detail the eigenvalues and eigenf unctions 
for conserving- and nonconserving-parity linear potentials. 

In the presence of a time-independent pseudoscalar potential the 1+1 
dimensional time-independent Dirac equation for a fermion of rest mass m 
reads 

(cap + (3mc^ + i; = Eip (1) 

where E is the energy of the fermion, c is the velocity of light and p is the 
momentum operator, a and (3 are Hermitian square matrices satisfying the 
relations = (3"^ = 1, {a, f3} = 0. From the last two relations it follows that 
both a and /3 are traceless and have eigenvalues equal to ±1, so that one 
can conclude that a and (3 are even-dimensional matrices. One can choose 
the 2x2 Pauli matrices satisfying the same algebra as a and f3, resulting in a 
2-component spinor ip. The positive definite function l^/^p = ip'^i/j, satisfying 
a continuity equation, is interpreted as a probability position density and its 
norm is a constant of motion. This interpretation is completely satisfactory 
for single-particle states [2H]- It is worth to note that the Dirac equation is 
covariant under x —x if V{x) changes sign. This is because the parity 
operator P = exp{ie)Poa3, where e is a constant phase and Pq changes x 
into —X, changes sign of a and P'j^ but not of p. Using a = ai and P = a^, 
P'y^ = (72 and provided that the spinor is written in terms of the upper and 
the lower components 

the Dirac equation decomposes into : 

— (^E — mc^^ ipj^ = ihcil)'_ + iVip- 

(3) 

— (^E + mc^^ ifj- = ihctjj'j^ — iVifjj^ 

where the prime denotes differentiation with respect to x. In terms of ip+ and 
?/'_ the spinor is normalized as + IV'-P) = so that ip+ and ?/'_ 

are square integrable functions. It is clear from the pair of coupled first-order 
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differential equations that ■?/'+ and have definite and opposite parities 
if the Dirac equation is covariant under x — > —x. In the nonrelativistic 
approximation (potential energies small compared to mc^ and E ^ mc^) Eq. 
Q loses all the matrix structure and becomes 

where V^^j = V'^ I2m(? + h/2mcV'. Eq. (j3]) shows that if of order 
v/c « 1 relative to ■?/'+ and Eq. © shows that ip+ obeys the Schrodinger 
equation with the effective potential K^j. It is noticeable that this peculiar 
(pseudoscalar-) coupling results in the Schrodinger equation with an effective 
potential in the nonrelativistic limit, and not with the original potential itself. 
The form in which the original potential appears in the effective potential, 
the V'^ term, allows us to infer that even a potential unbounded from below 
could be a confining potential. This phenomenon is inconceivable if one starts 
with the original potential in the nonrelativistic equation. 

The coupling between the upper and the lower components of the Dirac 
spinor can be formally eliminated when Eq. is written as second-order 
differential equations: 



2m 

where 



i'l + Kif fa = (6) 



= ^^f^ (7) 

Vif = ^—±-^V' (8) 
2mc2 2mc ^ ' 



These last results show that the solution for this class of problem consists 
in searching for bounded solutions for two Schrodinger equations. It should 
not be forgotten, though, that the equations for ipj^ or are not indeed in- 
dependent because the effective eigenvalue, -Ee//, appears in both equations. 
Therefore, one has to search for bound-state solutions for V^jj and V~jj with 
a common eigenvalue. 
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Now let us consider a pseudoscalar potential in the form 



hcg _ , 

V = muocx H (9) 

X 

where uj and g are real parameters. Then the effective potential becomes the 
singular harmonic oscillator 



Vtff = Ax^ + ^ + C^ (10) 



where 



A = -mio^ 
2 

= ^(g^l) (11) 

C± = huj(^g±^ 

It is worthwhile to note at this point that the singularity at x = never 
menaces the fermion to collapse to the center jl] because in any condition 
B± is never less than the critical value Be = —h /8m. Furthermore, to must 
be different from zero, otherwise there would be an everywhere repulsive 
effective potential as long as the condition 5+ < and i?_ < is never 
satisfied simultaneously. Therefore, the parameters of the effective potential 
with UJ ^ fulfill the key conditions to furnish spectra purely discrete with 
infinite sequence of eigenvalues, regardless the signs of u and g. Note also 
that the parameters of the effective potential are related in such a manner 
that the change cu -co induces the change V^jr — *• V^j — 2C± {C± —C±) 
whereas under the change g —g one has V^jr V^jr — 2C^ {B± — > B^ 
and C± — C*=f)- '^^^ combined transformation uj —uo and g —g has 
as effect Vg^j — > K//, meaning that the effective potential for ■?/'+ transforms 
into the effective potential for ?/'_ and vice versa. In Figure 1 is plotted the 
effective potential for some illustrative values of g. From this Figure one can 
see that in the case > 1 both components of the Dirac spinor are subject to 
the singular harmonic oscillator, for |5f| = 1 one component is subject to the 
harmonic oscillator while the other one is subject to the singular harmonic 
oscillator, for \g\ < 1 one component is subject to a bottomless potential well 
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and the other one is subject to the singular harmonic oscillator, in the case 
g = both components are subject to harmonic oscillator potentials. 
Defining 



f = ^ (12) 



and using (jH))-® one obtains the equation 



Now the prime denotes differentiation with respect to C,- The normalizable 
asymptotic form of the solution as ^ ^ oo is e^^/^. As ^ — > 0, when the term 
dominates, the regular solution behaves as where s is a nonnegative 
solution of the algebraic equation 

s{s - 1) - 2mB±/h'^ = (14) 

viz. 

.4(l±^l + ^)>0 (15) 

If Bj- > there is just one possible value for s (that ones with the plus sign 
in front of the radical) and the same is true for B± = Be when s = 1/2, but 
for Be < B± < there are two possible values for s in the interval < s < 1. 
If the singular potential is absent (-B± = 0) then s = or s = 1. The solution 
for all ^ can be expressed as V'±(0 = ^''^^e~^/^w;(^), where w is solution of 
the confluent hypergeometric equation [211 

^w" + {b- ^)w' -aw = (16) 

with 
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■mB± 1 



1 2m 



i^eff - C±) 



^± = 



(13) 



h 1 /2m 



s + 1/2 



(17) 
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Then w is expressed as i-Fi(a,6, ^) and in order to furnish normahzable %l)±, 
the confluent hypergeometric function must be a polynomial. This demands 
that a = —n, where n is a nonnegative integer in such a way that b,^) is 

proportional to the associated Laguerre polynomial L^~^{^), a polynomial of 
degree n. This requirement, combined with the top line of (fT7|). also implies 
into quantized effective eigenvalues: 

Eeff = (2n + s + ^^h\uj\ + C±, n = 0,1,2,... (18) 
with eigenfunctions given by 

^±(x) = iv±r/v«/2Lr^/2(^) (19) 

On the other hand, s = and s = 1 for the case B± = and the associated 
Laguerre polynomial and are proportional to (v^) 

^~^/'^H2n+i (v^); respectively j^Hj- Therefore, the solution for the harmonic 
oscillator can be succinctly written in the customary form in terms of Hermite 
polynomials: 



Eeff = [n + -j h\iu\ + C±, n = 0,l,2,... (20) 

^±(x) = iV±e-«/2if„ (v^) (21) 

Note that the behaviour of iIj± at very small ^ implies into the Dirichlet 
boundary condition {ip±{0) = 0) for s 7^ 0. This boundary condition is 
essential whenever B± 7^ 0, nevertheless it also develops for B± = when s = 
1 but not for s = 0. Since V^j is invariant under reflection through the origin 
{x — >■ —x), eigenfunctions with well-defined parities can be found. However, 
the eigenfunctions expressed in terms of associated Laguerre polynomials, due 
to the presence of the term x*, are restricted to the half-line x > in the event 
of s is not an integer number. This drawback can be circumvented by taking 
X by Then, there still an even eigenfunction such as in the case when s 
is an even number. New eigenfunctions can be found by taking symmetric 
and antisymmetric linear combinations of the eigenfunctions expressed in 
terms of eigenfunctions defined on the positive side of the x-axis. These new 
eigenfunctions possess the same effective eigenvalue so that there is a two- 
fold degeneracy. This degeneracy in a one-dimensional quantum-mechanical 
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problem is due to the fact that the eigenfunctions, whether they are either 
even or odd, disappear at the origin. These questions do not occur when the 
eigenfunctions are expressed in terms of Hermite polynomials because they 
are well-behaved on the entire x-axis and have both even and odd parities 
from the beginning. 

The necessary conditions for confining fermions in the Dirac equation 
with the potential have been put forward. The formal analytical solutions 
have also been obtained. Now we move on to consider a survey for distinct 
cases in order to match the common effective eigenvalue. As we will see this 
survey leads to additional restrictions on the solutions, including constraints 
involving the nodal structure of the upper and lower components of the Dirac 
spinor. The effective potentials for a few specific cases are already plotted in 
Figure 1 and all the others can be obtained by making the transformation 
ijj — > —ijj and g — >■ —g, as mentioned before. 

Case K-\g\>\ 

This is the situation where B± > and for (7>1((7<— l)s=|(7| for 
^// (^e//)) s = l + \g\ for V^jr (K//)- The upper and lower components 
of the Dirac spinor share the same eigenvalue if the quantum numbers satisfy 
the relation 



where e{uj) and e{g) stand for the sign function and = n (n_) is related 
to ip+ (V'-)- The solutions are 



= n + 



ejuj) - e{g) 
2 



(22) 




(0 



(23) 




2n + l- 



with 



n > 



e{g)-e{u) 
2 



(24) 
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and the further proviso that n > if e^u) > s{g). Note that the subscripts 
±£{g) in T/^ is a sequel of the transformation V^jr V^j—2Czf under g — > —g. 
Note also that the constraint ()22p asserts that the number of nodes of ip^ 
and ip- just differ by ±1. 



Case B — \g\ = 1 

In this situation 5+ = = 0) and s = or 1 for V^j- (K^j), and 
> (5+ > 0) and s = 2 for V'fj (V^^jf) when ^ = +1 = -1). 
It follows from (jl8j) and (j2(jp that the quantum numbers satisfy n_^(^g) = 
|n+e(g) — [2 — e{g)e{uj)]^ /2. Note that n_^i^g) is an integer number if and 
only if n^s(_g-) is an odd number. Making the substitution n —>■ 2n + 1 one 
can get the same constraint as in Case A, Eq. (j^^ . and the solutions can be 
written as 



(25) 



E. 



eff 



2n + 2 



h\uj\ 



and the same remarks for the Case A involving ip±e{g), as well as the nodal 
structure of the eigenfunctions, apply in Case B. 



Case C - < 1^1 < 1 

Here the solutions break in two sets regarding the two distinct possibilities 
of s for B+ <0 (5_ < 0) for V^^jf (V-jf) when < ^ < 1 (-1 < ^ < 0), 
viz. s = \g\ OT 1 — \g\. On the other side, > (5+ > 0) for V^j- (K//) 
and s = 1 + \g\. For the first set, s = for V^j- (K//)) the constraint 
involving the quantum numbers as well as the solutions for the effective 
eigenvalues and eigenfunctions are the very same as those ones presented in 
the Case A. For the second set, s = 1 — for V^j (K//)' the quantum 
numbers satisfy ^-^(g) = 'n.e(g) — \g\ + e{g)e{uj)/2. This second set provides 
n equal to an integer number on the condition that \g\ = 1/2 and this very 
particular case is already built in the first set. 
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Case D - ^ = 

This is the situation where B± 
satisfy 



= and s = or 1. The quantum numbers 



ri- — n + e{u!) 



(26) 



and the solutions are 



(27) 



Eeff = n + 



2 



hloi) 



with 



n > —s{uj) 



(28) 



and the further proviso that n > if£:(a;) = +1. 

The preceding analyses shows that the effective eigenvalues are equally 
spaced with a step given by 2h\u!\ when at least one the effective potentials is 
singular at the origin. It is remarkable that the level stepping is independent 
of the sign and intensity of the parameter responsible for the singularity of 
the potential. Wc also note that there is a continuous transition from the 
Case A to the Case C, despite the appearance of the Hcrmitc polynomial in 
the Case B. It should be remembered, though, that Hcrmite polynomials can 
be seen as particular cases of associated Laguerre polynomials. When both 
effective potentials become nonsingular {g — 0) the step switches abruptly 
to h\u\. There is a clear phase transition when (7 — > due the disappearance 
of the singularity for both components of the Dirac spinor. In the limit as 
— > the Neumann boundary condition, in addition to the Dirichlet bound- 
ary condition always present for 51 7^ 0, comes to the scene. This occurrence 
permits the appearance of even Hermite polynomials and their related eigen- 
values, which intercalate among the pre-existent eigenvalues related to odd 
Hermite polynomials. The appearance of even Hermite polynomials makes 
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ipi^) 7^ and this boundary condition is never permitted when the singular 
potential is present, even though g can be small. You might also understand 
the lack of such a smooth transition by starting from a nonsingular potential 
{g = 0), when the solution of the problem involves even and odd Hermite 
polynomials, and then turning on the singular potential as a perturbation of 
the g = potential. Now the "perturbative singular potential" by nature 
demands, if is either attractive or repulsive, that ip{0) =0 so that it naturally 
kills the solution involving even Hermite polynomials. Furthermore, there is 
no degeneracy in the spectrum for g = 0. 

In all the circumstances, at least one of the effective potentials has a well 
structure and the highest well governs the value of the zero-point energy. 
Matching the formal solutions for a common effective eigenvalue has imposed 
additional restrictions on the allowed eigenvalues and relations between the 
number of nodes of the upper and lower components of the Dirac spinor have 
been obtained. A sharp limitation occurred in the Case B = 1) when all 
the solutions involving even Hermite polynomials have been suppressed. In 
the Case D ((7 = 0), a case which cannot be obtained as a limit of the preced- 
ing ones, the solutions of the generalized Dirac oscillator were obtained. 
The Dirac eigenvalues are obtained by inserting the effective eigenvalues in 
((Tj). One should realize that the Dirac energy levels are symmetrical about 
E = 0. It means that the potential couples to the positive-energy component 
of the spinor in the same way it couples to the negative-energy component. 
In other words, this sort of potential couples to the mass of the fermion 
instead of its charge so that there is no atmosphere for the spontaneous pro- 
duction of particle-antiparticle pairs. No matter the intensity and sign of the 
coupling parameters, the positive- and the negative-energy solutions never 
meet. There is always an energy gap greater or equal to 2mc^, thus there is 
no room for transitions from positive- to negative-energy solutions. This all 
means that Klein's paradox never comes to the scenario. 

Figures 2-5 illustrate the behaviour of lip-l"^ and lipl"^ = \ip^\'^+\ip-\'^ 

for the positive-energy solutions of the ground-state contemplating all the 
distinct classes of effective potentials plotted in Figure 1. The relative nor- 
malization of ip+ and ip- is obtained by substituting the solutions directly 
into the original first-order coupled equations (jH)). Comparison of these Fig- 
ures shows that {ip^l is larger than (for E > 0) and that the fermion 
tends to avoid the origin more and more as l^^l increases. A numerical calcu- 
lation of the uncertainty in the position (with 'm = u! = c = h = l) furnishes 
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0.844, 1.130, 1.346 and 1.528 for g equal to 0, 1/2, 1 and 3/2, respectively. 

In conclusion, we have succeed in searching for Dirac bounded solutions 
for the pseudoscalar potential V = mujcx + Ticg/x. The satisfactory com- 
pletion of this task has been possible because the methodology of effective 
potentials has transmuted the question into Sturm-Liouville problems with 
effective quadratic plus inversely quadratic potentials for both components of 
the Dirac spinor. As stated in the first paragraph of this work, the anomalous 
magnetic-like interaction in the four- dimensional world turns into a pseu- 
doscalar interaction in the two-dimensional world. The anomalous magnetic 
interaction has the form —i^(3a.'V4>{r), where /i is the anomalous magnetic 
moment in units of the Bohr magneton and is the electric potential, i.e., 
the time component of a vector potential Therefore, besides its impor- 
tance as new analytical solutions of a fundamental equation in physics, the 
solutions obtained in this paper might be of relevance to the confinement of 
neutral fermions in a four- dimensional world. 
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Figure captions 

Figure 1 - Effective potentials for V — mwcx + Ucg/x with a; > 0, 51 > 0. 
The solid lines are for K/j, the dashed lines are for V^//- a) = 3/2 ; b) 
g — 1 ] c) g — 1/2 ; d) g — {m — cu — c — h = 1). 

Figure 2 - (full thin hue), (dashed line), I'j/'+P + K^^'-P (full thick 

line), corresponding to positive-gronnd-staXe energy for the potential V — 
mujcx + hcg/x with g = 3/2 and m = uj = c = h= l. 

Figure 3 - The same as in Figure 2 with g = 1. 

Figure 4 - The same as in Figure 2 with g — 1/2. 

Figure 5 - The same as in Figure 2 with g — 0. 
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